Implicitly Restarted Arnoldi Methods and Eigenvalues of the Discretized Navier Stokes Equations
نویسندگان
چکیده
Implicitly restarted Arnoldi methods and eigenvalues of the discretized Navier Stokes equations. Abstract We are concerned with nding a few eigenvalues of the large sparse nonsymmetric generalized eigenvalue problem Ax = Bx that arises in stability studies of incompressible uid ow. The matrices have a block structure that is typical of mixed nite-element discretizations for such problems. We examine the use of shift-invert and Cayley transformations in conjunction with the implicitly restarted Arnoldi method along with using a semi-inner product induced by B and puriication techniques. Numerical results are presented for some model problems arising from the ENTWIFE nite-element package. Our conclusion is that, with careful implementation, implicitly restarted Arnoldi methods are reliable for linear stability analysis.
منابع مشابه
Implicitly restarted Arnoldi with puri cation for the shift-invert transformation
The need to determine a few eigenvalues of a large sparse generalised eigenvalue problem Ax = Bx with semi-positive deenite B arises in many physical situations, for example, in a stability analysis of the discretised Navier-Stokes equation. A common technique is to apply Arnoldi's method to the shift-invert transformation, but this can suuer from numerical instabilities as is illustrated by a ...
متن کاملImplicitly restarted Arnoldi with purification for the shift-invert transformation
The need to determine a few eigenvalues of a large sparse generalised eigenvalue problem Ax = λBx with positive semidefinite B arises in many physical situations, for example, in a stability analysis of the discretised Navier-Stokes equation. A common technique is to apply Arnoldi’s method to the shift-invert transformation, but this can suffer from numerical instabilities as is illustrated by ...
متن کاملKrylov Methods for the Incompressible Navier-Stokes Equations
Methods are presented for t ime evolution, steady-state solving and linear stability analysis for the incompressible Navier-Stokes equations at low to moderate Reynolds numbers. The methods use Krylov subspaces constructed by the Arnoldi process from actions of the explicit Navier-Stokes right-hand side and of its Jacobian, wi thout inversion of the viscous operator. Time evolut ion is performe...
متن کاملSome new restart vectors for explicitly restarted Arnoldi method
The explicitly restarted Arnoldi method (ERAM) can be used to find some eigenvalues of large and sparse matrices. However, it has been shown that even this method may fail to converge. In this paper, we present two new methods to accelerate the convergence of ERAM algorithm. In these methods, we apply two strategies for the updated initial vector in each restart cycles. The implementation of th...
متن کاملA Harmonic Restarted Arnoldi Algorithm for Calculating Eigenvalues and Determining Multiplicity
A restarted Arnoldi algorithm is given that computes eigenvalues and eigenvectors. It is related to implicitly restarted Arnoldi, but has a simpler restarting approach. Harmonic and regular RayleighRitz versions are possible. For multiple eigenvalues, an approach is proposed that first computes eigenvalues with the new harmonic restarted Arnoldi algorithm, then uses random restarts to determine...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 1997